Evidence for Two Time Scales in Long SNS Junctions 
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We use microwave excitation to elucidate the dynamics of long superconductor / normal metal / 
superconductor Josephson junctions. By varying the excitation frequency in the range 10 MHz - 40 
GHz, we observe that the critical and retrapping currents, deduced from the dc voltage vs. dc current 
characteristics of the junction, are set by two different time scales. The critical current increases 
when the ac frequency is larger than the inverse diffusion time in the normal metal, whereas the 
retrapping current is strongly modified when the excitation frequency is above the electron-phonon 
rate in the normal metal. Therefore the critical and retrapping currents are associated with elastic 
and inelastic scattering, respectively. 

PACS numbers: 74.45.+c,74.40.+k,73.23.-b,72.30.+q 



Despite the large amount of work devoted to supercon- 
ductor / normal metal / superconductor (SNS) Joseph- 
son junctions, starting from the pioneering work of J. 
Clarke their dynamics remains an open issue. The 
main raison is that the time dependence of the phase 
difference between the two superconductors, well under- 
stood in tunnel junctions, interferes here with that of the 
out-of-equilibrium quasi-particles in the normal metal. 
For instance, the origin of the retrapping current 7^, the 
current below which the junction switches back from the 
dissipative to the non dissipative state, is still under de- 
bate 

One way to probe the junction dynamics is to mea- 
sure its dc response to an ac excitation. By varying 
the frequency of the excitation, the relevant time scales 
can be revealed. A BCS superconductor driven out- 
of-equilibrium by microwave radiation experiences in- 
creased superconductivity. This counterintuitive phe- 
nomenon is known as Dayem-Wyatt effect [4^, , first 
observed in microbridges, and then in thin films and 
strips 0, 0, H, S It is based on the fact that the 
BCS gap A is related in a self-consistent way to the 
energy distribution of the quasiparticles /(e) by: 1 = 
A °°p(e)[l — 2f{e)\de where A measures the strength 
of the electron-phonon interaction and p{e) is the BCS 
density of states, which depends on A. The thermody- 
namic gap A(T) corresponds to the equilibrium Fermi- 
Dirac statistics at temperature T; driving the system 
out-of-equilibrium with an ac current modifies the en- 
ergy distribution of the excitations, and thus the gap. 
To put the quasiparticles significantly out of equilibrium 
clearly demands enough power, but also enough speed: 
inelastic processes tend to restore equilibrium, and as a 
consequence enhanced gap is observed only when the ex- 
citation frequency / is greater than the inverse of the 
dominant inelastic process, electron-phonon in usual su- 
percoirductors. From the above equation, A(T) is found 
to be fairly insensitive to T at low temperature whereas it 
varies strongly close to the critical temperature Tc- This 



reflects the sensitivity of A(r) to a small change in the 
distribution function, and thus it translates into the fact 
that microwave pumping is effective only very close to 
Tc, experimentally for T within of few percent of Tc 0]. 

Hybrid structures made of superconductors and nor- 
mal metals offer a rich extension of the physics of su- 
perconductivity. In the SNS geometry, due to coherent 
Andreev reflections at the N/S interface, a finite gap A, 
the so-called mini-gap, develops inside the normal metal. 
In a way similar to the BCS gap equation, A is deter- 
mined by the Andreev bound states as well as by the 
distribution function of the quasiparticles. It is thus nat- 
ural to expect an enhanced induced superconductivity 
by ac irradiation 'l2'|. However, two facts make the 
physics of SNS samples more complex than that of su- 
perconductors: i) the distribution functions involve not 
only the inelastic time but also the diffusion time of 
quasiparticles along the N part, since Andreev pairs dif- 
fuse from the N/S interfaces to the center of the N part; 
ii) the density of states is affected by the rf field through 
the phase of the Andreev pairs being modified by the ac 
vector potential. 

We report an experimental investigation of the influ- 
ence of a microwave excitation on both the retrapping 
current Ir and the critical current /c, the maximum su- 
percurrent that can flow through the junction. 

We observe that a microwave of high enough frequency 
f > fc enhances Ic- We relate /c, temperature indepen- 
dent, to the inverse diffusion time through the normal 
wire. On the other hand, the switching back to the su- 
perconducting state is modified when f > fr, where fr is 
of the order of the electron-phonon scattering rate tJI^^ 
in the normal metal, and has the same, strong temper- 
ature dependence fr — T^lp^ oc {fr <C fc in our 
temperature range). 

The samples - The samples were fabricated by double 
angle deposition through a suspended mask, formed by a 
bilayer of 400 nm PES and 60 nm Si3N4. After drawing 
the mask by e-beam lithography, we etched with a SFg 
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FIG. 1: Left: Schematics of the experimental setup; Right; 
SEM picture of sample S. 



RIE the Si3N4 layer and then the PES with a O2 plasma. 
The undercut is about 0.5 /im [13]. The junctions were 
made by depositing a first layer of Al (9 nm), here the 
normal metal, and a second, overlapping layer of Nb (50 
nm), the superconductor. The deposition was made at 
a base pressure of 10~^ mbar. The critical temperature 
of the Nb film ranges from 6.8 to 8 K depending on the 
evaporation conditions. 

We have measured six samples of two different lengths: 
short samples (simply called 'S' below) with length L = 
330 nm, width w — 300 nm and normal state resistance 
Rn ~ 6 fl, and long samples ('L'), with L = 780 nm, 
w = 320 nm and Rn = 14.5 ft. We deduce the diffusion 
coefficient D ~50cm^s~^, and the diffusion times td — 
L^/D, with Tjj^ ~ 8.2 GHz for the long samples and 46 
GHz for the short samples. The different samples of the 
same length gave similar results. All the measurement 
were performed above 1.4 K when the Al is normal. We 
observe no evidence for superconducting fluctuations in 
Al above this temperature [l3| . 

Experimental setup - The experimental setup is 
sketched in Fig. [TJ The samples are dc current-biased 
through a 100 kfi resistor connected, after filtering, to 
the dc port of a bias tee, whereas a microwave genera- 
tor is connected to the rf port. The bias tee is essential 
for this experiment because it provides a wideband, good 
coupling of the sample to the rf source. The sample is 
placed at the end of a 50 stainless steel coaxial line 
and sits in the liquid of a pumped He* cryostat. The 
temperature is varied by adjusting the bath pressure. 
When it is superconducting, the sample has a very small 
impedance, thus the rf generator of output impedance 
Rq ~ 50ri acts as a current source. This is not strongly 
modified when the sample is in the normal state since 
its resistance is ~ 0.1 — 0.3i?o- At high frequency the 
sample has a complex impedance that can be of order 
i?.o 15]. However, it is wire bonded to the microstrip on 
the sample holder. We have measured the inductance in 
series with the sample, L = 3.3 nH, by refiectometry in 
the range 0.1-5 GHz. This is the dominant impedance 
at high frequency (Lw ^ 200 O at 10 GHz), so the sam- 
ple is still current-biased at high frequency, but the ac 
current is reduced by a factor Rq/{Lu}) that we take into 
account, together with the attenuation of the cables that 
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FIG. 2: (color online) dc voltage vs. dc current measured on 
sample S with a 10 MHz (a) or 100 MHz (b) rf excitation 
at several powers : from right to left, P=-66 dBm, -52dBm, 
-48dBm and -46dBm. The dotted line corresponds to no rf 
excitation. The temperature is T = 1.5 K. 



we have measured separately, to estimate the rf power at 
the sample level. We observe several regimes depending 
on the frequency of the microwave: at low frequency the 
voltage follows adiabatically the dc measured V{I) curve; 
at intermediate frequency, the V{I) is modified and the 
measured critical current is lowered by the presence of 
the ac excitation; at high enough frequency we observe a 
strong enhancement of the critical current. 

The low and intermediate frequency regimes - A dc 
V{I) characteristics for sample S at T=1.5K with no rf 
excitation is shown on Fig. [21 dotted line. One observes 
a hysteresis (for T < 2.23 K, see FigIS]) from which we 
deduce the critical current J° and the retrapping current 

with no excitation. Let us first consider what happens 
when a slow ac current of amplitude lac is added to the 
dc current. Increasing the dc current I from zero we ex- 
pect the sample to jump to the N state when I + Iac = Ic 
(for small enough lac)- Decreasing / from the N state, 
the sample stays normal as long as I~Iac > ll-- Thus for 
lac < Iq ~ J cxpccts One jump at Ic = Ic — lac, i-e. 
the apparent critical current Ic decreases with increased 
rf power. When lac > Ic^^r, the sample will cycle from 
S to N and vice- versa, as long as 1° — lac < I < Ir + lac- 
Thus in this regime the V(I) characteristic should have 
two steps at I^ = Ic - lac and I^ = 1° + lac- The 
two regimes with one step for small lac and then a dou- 
ble step for larger lac is exactly what we observe at 10 
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FIG. 3: (color online) /c(r) (+: rf off, Q: rf on) and Ir{T) 
(x: rf off, A: rf on) with and without rf excitation (39 GHz, 
-17.4dBm) on sample S; dotted line: theory for 9.4 x Ic{T). 
Inset: symbols: fr vs. T for sample S; solid line: electron- 
phonon scattering rate for Al vs. T. 

MHz, see Figl2ja), the position of the steps varying as 
discussed above. Thus the sample follows adiabatically 
the 10 MHz excitation. When we increase the excita- 
tion frequency, the double steps disappear, see Fig[21[b). 
More precisely, the measured critical current still verifies 
Ic — Ic ~ lac even at high rf power. Therefore when the 
sample becomes normal, it never jumps back into the su- 
perconducting state even though the current I{t) spends 
some time below the dc retrapping current This dy- 
namical effect can be understood if we suppose that the 
hysteresis in the V{I) curve is at least partially due to 
heating of the electrons by the dissipative current, as sug- 
gested in : the sample can return to the N state only if 
its temperature is low enough, i.e. when there is not too 
much Joule power dissipated by the current, which at low 
frequency occurs at a current / = /° < 7°. In our case, 
if the frequency is larger than the inverse thermalization 
time of the electrons, the instantaneous current may be 
below 7° but electrons are still hot and the sample is in 
the normal state. The thermalization of the electrons 
occurs via phonon emission. In the inset of FigO we 
present the temperature dependence of the frequency fr 
corresponding to the transition from two-step to one-step 
I(V) curves, as well as that of the electron-phonon scat- 
tering rate for Al, rjl'p^ — AT'^, with A being 2.1 times 
larger than the theoretical value, in good agreement with 
previous measuremenst [l^. This strongly suggests that 
the retrapping is related to energy relaxation by electron- 
phonon interaction. 

The high frequency regime - We will now consider the 
high frequency regime. The V{I) characteristics has only 
single jumps. To measure Ic we current-bias the junction 
with an 83 Hz linear ramp and measure the time spent 
before dissipation. In Fig. [3] we show 7c(T) and Ir{T) 
for sample S with and without ac excitation (39 GHz, 



-17.4dBm). We observe that Ic is strongly increased by 
the rf excitation over the whole temperature range. This 
is similar to what has been observed in Pb/Cu/Pb junc- 
tions [13] but in strong contrast with what is observed 
in superconducting strips, where the enhanced supercon- 
ductivity is restricted to a very narrow region around 
Tc, of order of a few mK for Al. The dotted curve of 
Fig. [3] corresponds to ^ 9.4 times the theoretical ex- 
pectation for 7c (T) in a long SNS junction [l3| (there 
are no adjustable parameters). We think that the factor 
9.4 is due to the non zero superconducting coupling con- 
stant in the Al. We will now focus our discussion on the 
critical current and discuss the retrapping current later. 
We report in Fig. |4] the effect of the rf current on Ic 
for various frequencies. We observe that below a certain 
frequency fc the critical current is decreased by the rf 
excitation, whereas for / > fc, Ic is increased. We ob- 
serve an enhancement of Ic by almost a factor two on L 
samples irradiated at 39 GHz. Ic varies continuously vs. 
rf power (monotonically up to 33 GHz) up to a threshold 
at which the sample becomes suddenly normal. This be- 
havior is strikingly different from what has been reported 
in Sn/Au/Sn junctions and attributed to phase effects, 
where the increased critical current smoothly oscillates 
as a Bessel function of the rf voltage [lH. Our S and 
L samples show a similar behavior of Ic vs. rf power, 
but with a frequency fc that is larger for short samples: 
fc ^ 7 GHz - 0.9r^^ for L samples and fc ~ 17 GHz 
~ OAtjj^ for S samples. In ref. J^] is mentioned that the 
relevant time scale may be the "effective time-dependent 
Ginzburg-Landau relaxation time" t* proportional to the 
diffusion time: t* = (7r/2)^r£). 

In order to elucidate the mechanism that sets the value 
of fc, we have: i) varied the temperature, and ii) applied 
a magnetic field 77 perpendicular to the substrate. We 
observe: i) no noticeable effect of the temperature on fc, 
and ii) L and S samples show a non-monotonic behavior 
of fc{H), the field dependence being very temperature 
dependent, see inset of FigHl This suggests a competi- 
tion between the field dependence of the mini-gap and 
the magnetic field induced broadening in the quasiparti- 
cle density of states. 

The retrapping current - We now discuss the effect of 
the rf excitation on the retrapping current. We show 
in Fig. [3] the temperature dependence of 7,. with and 
without rf excitation. One clearly distinguishes two tem- 
perature ranges: at high temperature, T > T* ^ 2.23 
K, there is no hysteresis in the V{I) characteristic, so 
Ir = Ic is enhanced by the rf excitation. For T < T* , 
the value of Ir saturates and is barely modified by the rf 
excitation. Thus, as soon as Ir is a relevant quantity (i.e. 
when V{I) is hysteretic), its response to the rf excitation 
is radically different from that of Ic- We obtain similar 
results for samples S with T* ^ 1.6 K. 

The hysteretic behavior of SIS (I=insulator) Josephson 
junctions comes from the phase dynamics of the junction 
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FIG. 4: (color online) Critical current vs. rf power for dif- 
ferent frequencies: from bottom to top f=3, 6, 12, 21, 33, 36, 39 
GHz. Inset: symbols: renormalized cross-over frequency 
fc{H)/fc{0) vs. renormalized magnetic field H/Ho, where Ho 
corresponds to one flux quantum in the N part. Full squares: 
sample L at T = 1.6 K; full dots: sample S at T = 1.4 K; 
circles: sample S at T = 2 K. 



tivity in Ss'S occurs very close to Tc. In contrast, the 
effect of rf on a Josephson junction arises through the 
time dependence of the Josephson phase, which leads 
e.g. to Shapiro steps[2l|. Our situation is intermediate 
and a quantitative analysis requires a careful theoretical 
study of the interplay between out-of-equilibrium super- 
conductivity and phase dynamics in these junctions [22l |. 
However, using ac irradiation we performed a photon- 
assisted like experiment where we have identified two dif- 
ferent timescales for the critical and retrapping currents. 
Specifically, we showed that the critical current and the 
retrapping current are associated with elastic and inelas- 
tic scattering, respectively. 

We are grateful to J.-C. Cuevas for showing us some 
of his results prior to publication. We thank H. Bouchiat, 
M. Ferrier, J. Gabelli, S. Gueron, I. Petkovic for fruitful 
discussions and I. Petkovic for technical assistance in the 
sample preparation. This work was supported by ANR- 
08-BLAN-0087-01. 



that is highly sensitive to the electromagnetic environ- 
ment. The phase dynamics is usually described by the 
RSCJ model [iSj. Following this model and substituting 
the characteristic (vo ltag e relaxation) time, RC, by the 



diffusion time to 3, 



20[, one expects Ir cc \JTc^ in con- 
tradiction with our measurement. Another source of hys- 
teresis in SNS junction is of thermal origin: when in the 
normal state the sample dissipates Joule power R^I^- 
Because of the superconducting contacts, the electrons 
can relax their energy only by phonon emission whose 
cooling power is Pe-ph oc — T^, with the elec- 
tron temperature and T the phonon bath temperature. 
To switch back in the superconducting state, the cur- 
rent and thus the electron temperature, have to be such 
that / < /c(Te). At high temperature the cooling is very 
efficient and Ir — Ic- At low temperature the cooling 
is poor, with Tg almost independent of T, and Ir satu- 
rates to a value significantly lower than Ic- This qualita- 
tively reproduces our measurement. However, from our 
measurement of Ic{T) we find that Ir should saturate at 
~ 7 /iA (i.e. Tg — 2.9 K). This value is 2.5 times smaller 
than what we measure. More importantly, Ir should be 
enhanced by the rf excitation as Ic is, and it is not. 

Conclusion - SNS systems are intermediate between 
weakened (constricted) superconductors (Ss'S geometry, 
s'=weak superconductor) and Josephson junctions with 
an insulating barrier, and one is tempted to interpret 
their behavior in terms of these two limits. In constricted 
superconductors one is dealing with usual bulk supercon- 
ductivity except that the populations arc driven out of 
equilibrium by the rf pumping. Enhanced superconduc- 
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